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Abstract. - We investigate dissociation reactions of loosely bound and spatially extended 
three-body systems. We formulate a practical method for simultaneous treatment of long-range 
Coulomb and short-range nuclear interactions. We use 6 He (n+n+a) and n Li (n+n+ 9 Li) as 
examples and study the two-neutron separation cross sections as functions of target and beam 
energy. Individual Coulomb and nuclear as well as interference contributions are also extracted. 



Introduction. - Nuclear halos were discovered by their surprisingly large reaction cross 
sections with ordinary target nuclei Jl|-^| . Most of the information about these unusual struc- 
tures are obtained by detailed fragmentation reaction studies [|s|-^0|. The understanding of 
these halo systems is based on an effective clustering into few-body structures reacting with 
the target. The main structure is essentially agreed upon as two or three-body systems with 
11 Be (n+ 10 Be) as the prototype of a two-body halo and 6 He (n+n+a) and n Li (n+n+ 9 Li) 
as three-body halos Q . 

The reaction description is much less advanced and are mostly available as fragmentary 
computations where one quantity is investigated in one particular model with one set of 
parameters fi|, ^2| |l6| . Systematic and consistent calculations of many observables within 
one model are highly desired but scarce 0. The difficulties in an accurate treatment arc 
substantial, since even when the constituent clusters and the target are oversimplifyingly 
considered to be point-like, the reactions still involve three and four-body systems for two and 
three-body halos, respectively. Approximations are therefore inevitable. 

One of the major problems in obtaining fragmentation cross sections is to incorporate 
Coulomb and nuclear interactions in the same numerical procedure J5,18|. This has been 
attempted for two-body projectiles [^9[^0| while three-body projectiles at best are treated as 
effective two-body systems JpLflq} or by computing the Coulomb and nuclear contributions in 
independent models |5|, p^, |lS|, |2^| . This situation is rather unsatisfactory, since breakup 
reactions of nuclear halos are dominated by nuclear and Coulomb interactions for light and 
heavy nuclear targets, respectively p, 18,p2[. 

Significant improvements can be expected to be difficult, especially for three-body pro- 
jectiles, since a mixture of long-range and short-range interactions is involved. The efforts 
may also be very rewarding first by allowing the necessary treatment of specific reactions and 
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second due to the general nature of the problem and the derived interest from other subfields 
of physics. We formulate in this report a practical method to include Coulomb and nuclear 
interactions simultaneously in investigations of breakup reactions of three-body halo nuclei. 

Model assumptions. - We consider a three-body halo system colliding with a target. We 
assume first that the intrinsic motion of the halo is slow compared to the relative projectile- 
target motion. The three-body system does not have time to adjust to the external field 
during the collision and the sudden approximation therefore applies. For point-like particles 
each halo particle, called the participant, can then interact with the target without disturbing 
the motion of the other two, called spectators. We can then treat the reaction as independent 
collisions and the total cross section is the sum of contributions from the three participants. 
Apart from energy and momentum conservation and overlaps between initial and final states 
of the spectators we are then left with three independent two-body problems. The corre- 
sponding interactions should then describe these two-body collisions to the required level of 
accuracy. We use the phenomenological optical model designed to describe elastic scattering 
and absorption from the elastic channel. 

The finite extension of both target and halo particles demands additional considerations, 
since simultaneous collisions of more than one halo particle then could be quite frequent in 
contradiction to our basic assumption. Therefore in addition to the detailed treatment of the 
essential contribution of the participant we use the simpler "black sphere" model to describe 
the smaller contributions from the spectators. If the spectators are able to pass without 
touching the target they are true spectators and otherwise they are counted as absorbed by 
the target and consequently removed from the final state. This is the optical model limit 
of very strong absorption for short-range potentials. This division into detailed treatment of 
participant and use of the black sphere for the spectators is only meaningful when the halo 
is larger than the combined sizes of the target and the participant. This is similar to the 
assumption used in the formulation for a weakly bound projectile ]23| . A better treatment of 
the collision between three-body halos and a target almost inevitably has to deal with more 
than three-body configurations in the final state or include properties of the intrinsic structure 
of the halo particles and the target. 

Method. - The finite extension of the projectile particles and the target destroy the clear 
division into participant and spectators, since the spectators may hit the target in the same 
collision where the participant is scattered or absorbed. We start by treating the interaction 
between the target and the projectile constituents in the black sphere model where the particle 
is absorbed inside a cylinder with the axis along the beam direction and left untouched outside 
this cylinder of a radius approximately equal to the target plus spectator radius. If P c , P n \ and 
P n 2 are the probabilities for the core, first and second neutron being inside this cylinder the 
reaction probability is the sum of the probabilities of finding all three constituents inside the 
cylinders {P n \Pn2Pc), plus the probability of finding two constituents inside (three terms like 
P%Pj(l — -Pfc))) phis the probability of finding only one constituent inside the cylinder (three 
terms like Pj(l — Pj)(l — Pk))- We then consider all possible combinations in the projectile- 
target interaction. The sum of all these terms give the reaction probability, that can be 
rewritten as P c + P n i(l — P c ) + fn2(l — ^c)(l — Pni)- Each term in this probability vanishes 
unless one particular projectile constituent is inside the cylinder. This is the constituent 
chosen as participant, while the other two are considered spectators. The optical model is 
then used for the participant-target interaction. 

Participant treatment. - We consider collisions for an initial velocity v, the corresponding 
momentum p, and total kinetic energy of EA, where A is the mass number of the projectile. 



E. Garrido, D.V. Fedorov and A.S. Jensen: Coulomb and nuclear breakup of three-body halo nuclei3 



The target (labeled 0) has charge Zq and mass mo, the participant (labeled i) has charge Zi 
and mass mj. We label the spectators by j and k, final state quantities by primes, relative 
two-body coordinates and momenta between particles i and k by and p ik and between i 
and the center of mass of j and k by r^jk and p i j k . The initial wave function is a product 
of the three-body wave function and a plane wave describing the relative halo-target motion. 
The final state wave function is a product of three terms, i.e. two distorted waves for the 
participant-target and the spectator-spectator motion and a plane wave for the relative motion 
of these two non- interacting two-body systems. The distorted waves are obtained by solving 
the Schrodinger equation with the appropriate two-body potentials. 

The cross section of the participant i has two contributions corresponding to elastic scatter- 
ing (diffraction) and absorption (stripping). The differential diffraction cross section reduces 
to a factorized form when the participant has spin or 1/2 and the target has spin |L7j 

dPo iJk dp' jk dp' 0i dp' 0i 
x (1 - meM^q-r^ k )\n 2 ) \M s (p 2jkl p' jk )\ 2 , (1) 



where p' Qi jfc = p i ^ k + Po( m i + m fc)/( m o + m i + m j + m k) in the rest frame of the projectile, 
^ is the initial three-body halo state, Sq — (p'ij k — Pijk)( m j + m k)/(rm + mj + m k ) is the 
momentum transfer into the relative participant-spectators motion and M s (p i j kl p'j k ) is the 
normalized overlap matrix element between initial and final state spectator wave functions. 

The first factor, d 3 a^ (p oi — > p' Qi )/dp' oi , is the differential cross section for the participant- 
target elastic scattering process. It is obtained by numerical computation of the phase shifts 
including simultaneously the nuclear and Coulomb potential. At small momentum trans- 
fers (large impact parameters) or large angular momenta the Rutherford cross section is 
approached as the effect of the short-range nuclear interaction then disappears. 

The second factor in Eq.(|l|) is constructed as one minus the probability for staying in 
the ground state after transfer of the momentum 5q in the reaction, i.e. we remove the 
probability for elastic scattering of the halo as a whole. This could alternatively be done by 
orthogonalizing the final state in the overlap matrix element to the initial three-body ground 
state. 

Collisions at sufficiently large impact parameters only produce virtual excitations corre- 
sponding to adiabatic motion [l8|]. The limiting impact parameter b a is determined by 
equating the reaction time 2b a /v with the time period 2nh/ B ps in the relative motion of the 
participant-spectators system with the corresponding binding energy B ps . The Coulomb in- 
teraction then transfers the momentum q a = ZoZie 2 p/ (b a AE) = ZyZie 2 /(ire) B ps /(hc) (7 + 
l)7~ 2 /3~ 2 , where j3 = v/c and 7 = — (3 2 . 

The energy transferred from target to participant, 5E = \/pq + — \l p 2 + m 2 ,, must 
be larger than the three-body separation energy B. When po and q = po Po' are parallel 
SE is maximized. For this geometry we find for small B compared to the target rest mass 
that SE = B implies that qc = q^c w By/ 1 + nri^c 2 jy\ which reduces to B/v in the non- 
relativistic limit. Thus q must be larger than qi, to produce dissociation, but on the other 
hand dissociation is not the necessary outcome for all q > qi- We exclude contributions from 
momentum transfer smaller than the largest of q^ and q a . 

The differential absorption (stripping) cross section, where the participant in the sense of 
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Fig. 1 - Two-neutron removal cross sections as functions of beam energy for fragmentation of 6 He 
(thick) and Li (thin) on a Pb-target. We show the total cross section (solid) as well as contributions 
from Coulomb (dashed), nuclear (dotted) and interference terms (dot-dashed). The curves for 6 He 
are lower than those of n Li. The experimental data are from flq-|j[fLo||, u Li (open circles) and 6 He 
(filled circles). 



the optical model is absorbed by the target, is obtained analogously |0 

*abs (P0i) \M s (p ujk ,p' jk )\ 2 , (2) 

where is the participant-target absorption cross section. The nine-dimensional differential 
cross section is now reduced to six, since the absorbed or stripped particle inherently is of 
no interest in the optical model description. The factorizations in both Eqs.(Q) and (||) are 
incomplete, since p oi via momentum conservation is related to p' oi - fe and p^jk- 

The total differential cross section for any reaction is now obtained by adding the two 
contributions from scattering and absorption to analogous contributions when particles j and 
k are the participants. All these contributions are given on an absolute scale and their relative 
weights are therefore determined in this model. We integrate over all unobserved momenta in 
Eqs.® and (§). 

Spectator treatment.. - The spectator-target interaction is treated in the black sphere 
model. This model can be implemented by considering the distance r ps between the partic- 
ipant and the spectator jL7|. When r ps is small participant and spectator both would for 
short-range interactions interact with the target, and therefore according to the black sphere 
model the spectator would be absorbed. The distance r ps is then determining if the spectator 
is either absorbed (distances smaller than r ps ) or undisturbed (distances larger than r ps ). The 
absorption distances r ps [jljj are related to the sizes of target and spectators and determined 
by f r ps = ( r2 )t + ( r2 )s + 2 fm 2 , where () is the measured mean square radius of target or 
spectator and 2 fm 2 is the square of the range of the nucleon interaction. These radii r ps 
can also be obtained from the parametrization in Eq.(28) of JTtJ with r = 1.26 fm, L30 fm 
and 1.45 fm for the targets Pb, Cu and C, respectively. In the same way when the charged 
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Fig. 2 - The same as Fig. |l| for a Carbon target. The experimental data are from 



participant is absorbed in the optical model sense, it must have been close to the target and 
close-lying spectators within the absorption distance must also be counted as absorbed. On 
the other hand when the charged participant is scattered small momentum transfer corre- 
sponds to a large impact parameter. These events only occur when the distance between 
participant and target is large and it is then very unlikely that spatially close-lying spectators 
are absorbed. All spectators are therefore counted as scattered. For large momentum transfer 
in the scattering process corresponding to impact parameters smaller than the sum of target 
and participant radii the spectator is analogously counted as absorbed. For two-neutron dis- 
sociation cross sections this division is irrelevant, since these contributions all are included in 
the cross section. 

Results.. We apply the procedure on the prominent nuclear three-body halos 6 He 
and n Li. The neutron-neutron and the neutron-core two-body interaction parameters are 
given in Jlj| |l7). The optical model parameters are from [^4| for the neutrons, from |25|] for 
a-particles and for 9 Li-particles from |25|] with range and diffuseness parameters fro m pc| . 
We furthermore drastically reduce the energy dependence of the real potential in |25fl to 
allow for the required huge beam energy variation, i.e. a? = —0.014. The measured core- 
target interaction cross sections are reproduced within error bars 0. The binding energy 
Bp S between the 9 Li and 4 He cores and the two neutrons must be introduced for the low 
momentum cutoff. We use the scaling relation in (27) to obtain B ps /B w 3 for 6 He and 1.4 
for n Li. The two-neutron dissociation cross sections are shown in fig. |l]as functions of beam 
energy for a Pb-target. The theoretical uncertainties are mostly due to the inaccurate optical 
model. Also the precise choice of cutoff parameter gives an uncertainty on the Coulomb 
part, especially important at higher energies. For Pb the adiabatic cutoff is the largest and 
therefore decisive, i.e. q a > qL- The n Li results are larger than those of 6 He because the 
Coulomb contribution roughly scales with the square of the projectile charge and the nuclear 
part increases with the projectile size. The small nuclear part has contributions from both 
neutron and core as participants. They both decrease up to about 200 MeV/nucleon then 
remaining roughly constant at higher energies. Discussions of the individual behavior of the 
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Fig. 3 - The same as fig. [l] for a Copper target. The experimental point for 11 Li is from 



many contributions are not possible here. The dominating Coulomb contributions decrease 
with energy first strongly and then much slower above 200 MeV/nucleon. The interference 
terms are negative and rather small, but significant at low energies, where the total cross 
section is below the Coulomb contribution. The measured values are reproduced at high 
energies but exceeded by a factor of two at low energy as already noticed in |l7j . Our results 
are consistent with estimates reported in . 

The trends for the individual terms are roughly the same for C as for Pb, see fig. g. The 
nuclear parts increase slowly with energy above 200 MeV/nucleon. The adiabatic cutoff is now 
the smallest and therefore is decisive, i.e. q a < q^- We observe a small increase of the total 
cross section with increasing energy. The nuclear contributions are completely dominating 
for a light target and the numerical values are therefore also very similar for both projectiles. 
Again we do not reproduce the (inconsistent) experimental data at lower energies. 

Coulomb and nuclear contributions are comparable for a medium heavy target, see fig. ^. 
The two cutoff parameters are now comparable and q a is decisive at small energies whereas 
qi, becomes larger at higher energies. The trends for the individual terms still remain the 
same as for C and Pb. Now the decreasing Coulomb contribution crosses the almost energy 
independent nuclear contribution at about 300 MeV/nucleon for 6 He and at 200 MeV/nucleon 
for n Li. The interference terms are still rather small except at low energies. 

Conclusion.. - We have formulated a method to compute dissociation cross sections of 
loosely bound three-body systems interacting with a mixture of short and long-range poten- 
tials. The uncertainties in the numerical results are dominated by the uncertainties in the 
two-body optical potentials for light targets and in the cutoff parameter for heavy targets. 
The already rather heavy computations seem to suffice as seen by comparing with the mea- 
surements. The present consistent and systematic calculations may prove useful as a guide 
for future experimental investigations. The dissociation cross sections are for example pre- 
dicted to increase as the beam energy decreases below 200 MeV/nucleon and slowly decrease 
or increase with increasing energy for higher energies, respectively for heavy and light targets. 



* * * 
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